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Surface states and the thermal Casimir interaction 
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Abstract - Using exact calculations, we elucidate the significance of the surface (bound) states 
for the thermal Casimir interactions for an Ising strip with a finite width. The surface state 
arises whenever an imaginary wavenumber mode appears in the spectrum of the transfer matrix, 
taken in the direction parallel to the edges of the strip. Depending on the boundary conditions, 
the imaginary modes emerge below or above the bulk critical temperature, or below the wetting 
temperature of a single surface with surface magnetic field. The bound states are responsible for 
the strong asymmetry of the Casimir forces between the super- and sub-critical regimes and for 
their sign. Our analysis uses the fact that the Casimir forces have two mathematical forms. We 
show that these very different representation are the same and in the process find the origin of 
the asymmetry. 
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It has been known for about half a century that the 
existence of quantum fluctuations in the electromagnetic 
vacuum is revealed by the existence of long-ranged, largely 
attractive, forces between objects pQ. For instance, in 
the geometry of parallel plates at a distance L apart, the 
wavenumbers in the direction perpendicular to the plates 
become discrete and this produces a change in the energy, 
which (per unit area) behaves as SE(L) cx hc/L 3 for suf- 
ficiently large distances, c is the speed of the light and h 
is a Planck's constant. Fisher and de Gennes [2] made 
the crucial observation that in condensed matter systems, 
1 say uniaxial classical ferromagnets and their analogues, 
the confining surfaces restrict order parameter fluctuations 
leading to to "entropic" forces between the surfaces. Scal- 
ing theory applied to these fluctuations predicts the exis- 
tence of forces per unit area in the critical region (at zero 
bulk field) of the form (in the following all free energies 
and forces are expressed in units of ksT) 



F C as=L- d Q{sgn{T-T c )L/£,) 



(1) 



where T c is the bulk critical temperature, £ is the bulk 
correlation length and the distance between the confining 
surfaces is L [2]. These forces are of significance on the 
micro- and nanoscale as shown in experiments on wetting 



films and colloidal systems [5HS]. 

The finite-size scaling function O is expected to depend 
on the geometry and on the boundary conditions (BCs) 
imposed by the confining surfaces. Mathematical consis- 
tency may allow small and large argument behavior in (fl} 
to be inferred, but to go further, detailed calculations of 
are needed. The main generally applicable techniques 
are Monte- Carlo simulations [7J and the de Gennes-Fisher- 
Upton local functional method [8] . The success of the first 
depends crucially on a thorough understanding of finite- 
size effects; the second is phenomenological in character, 
but with the benefit of containing no adjustable param- 
eters. In view of these matters, exactly solvable models 
have a special relevance, since they afford insights not oth- 
erwise available. 

In this Letter, we shall consider Casimir forces in the 
planar Ising model with a strip geometry. This encom- 
passes binary mixtures and lattice gases, with the surface 
magnetic fields playing the role of surface (differential) fu- 
gacities. Depending on the method of calculation used to 
obtain them, the results of the Casimir force come in two 
very different mathematical forms. One of them, a mode 
sum, has received little analytical attention so far. It has 
two advantages: firstly, the finite geometry imposes the 
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discretisation condition in this mode sum, as in the quan- 
tum Casimir case. The second advantage, and this is one 
of the results of this Letter, is that the strong asymme- 
try between the super- and sub-critical regimes receives 
for the first time a microscopic interpretation through the 
contribution of bound states, or surface modes. Their ex- 
istence in this context has been known for quite some time 
[9], but not the interpretation, which we provide in this 
Letter in terms of the Fisher-Privman finite-size scaling 
theory [12]. The other representation of the Casimir force 
for strips, this time as an integral, comes either from Au 
Yang and Fisher [T3E3] or from extension of the Schultz, 
Mattis and Lieb (SML) method QUO!], with the follow- 
ing advantages: it is ideal for taking the scaling limit, and 
also for investigating real analyticity in x of Q(x) from 
Further, although it is the obvious choice for numer- 
ical work, it gives no insight into the striking asymmetry 
mentioned above [T7] . The most recent contribution along 
this line explains how the boundary fields can be adjusted 
to produce both attractive and repulsive critical Casimir 
interactions [T7] . 

We now specify the model more precisely. As shown in 
Fig. [TJ we have an Ising strip with cyclic BCs in the (1, 0) 
direction. The energy for a configuration {er} of spins is 

M N 

e n ,m(w}) = -J2J2 Ki 

(2) 

m—l n—1 

M N-l 

^ y ^ y ^2^m,fi^m,?)41 ~t~ ^boundary- 



The term Ei, oun dary just concerns the rows n = 1 and 
n = N. This situation can be handled by a transfer matrix 
(TM) building the lattice in the (0, 1) direction, where the 
matrix has cyclic symmetry in the (1,0) direction; the 
SML technique [TS] may be applied. The alternative is to 
transfer in the (1,0) direction; cyclic symmetry no longer 
obtains, but the TM spectrum can still be found [9 , 18 - 20 
for the BCs which we consider. 

The free energy per unit length in the (1,0) direction, 
as derived from the transfer in the (0, 1) direction is |17) 



f(N) 



Air 



duj In 
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B x B 2 e 



(3) 



up to an additional term (N/2) ln(sinh 2K 2 ). Here, j(ui) 
is the Onsager function given by |21) 
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Fig. 1: (color online) An Ising strip with M columns and N 
rows. We are interested in the limit M — > oo with N < oo. 
Surface fields hi and hi are applied at the bottom and top 
edges; in the present paper we consider free boundaries hi = 
ha = and the (+, +) boundary conditions hi = ha = oo. 
In the model (a) the transfer matrix is along the (1,0) axis 
whereas in the model (b) along the (0, 1) direction. 



at the value +1, we have A\ = A 2 = e^cos((5*/2) and 
B x = B 2 = e-^sin(<5*/2). The angles 5* and 6' are ele- 
ments of the Onsager hyperbolic triangle [5T] , given by the 
formula sin 8' (uj) = sinh2i^i sin(w)/ sinh7(oj) and cos 8' = 
(cosh 2KI cosh 7 — cosh 2K\)/ (sinh 2isT| sinh 7) . The an- 
gle 5 (uj) has K\ and K 2 interchanged. Extracting from 
© the bulk contribution N f x = -N(l/4%) i(w)du>, 
which is easy, and an additional, N independent, surface 
contribution f s = —(I /Air) log [AiA 2 ] du 
gives the Casimir free energy: 



cosh 7 = cosh 2K X cosh 2K^-siah2K x sinh 2if| cos w, T Cas = — - / dtu In f 1 + e - 2N ^ ] tan 2 f 8{u)/i) ) 



(4) 

where the non-negative real branch is taken for real 
uj, K* is the dual coupling given by the involution 
8wh2KjSmh2K? = l,j = 1,2; Aj(uj) and B i (w) de- 
pend on the BCs. For free boundaries, we have A\ = 
A 2 = e* 1 cos((5'/2) and B 1 = B 2 = sm{8' /2), for (+, +) 
BCs with all spins on the bottom and top boundary fixed 



(5) 

where 8(ui) — Si(cu) + S 2 (lo), and Si(uj),i = 1,2 depend on 
the BCs. For free boundaries S(uj) = 8*(u>), whereas for 
(+,+) BCs S(uj) — 8'(u)). 

On the other hand, the transfer matrix V for 
the (1,0) direction has the diagonal form [SJ |T8 20 
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V = exp 



!£f=i7MA0)(2Xj^-l) 



The Xi 



are Fermion operators with anticommutation relations: 



Thus 2XjXj - 1 has 



[x v x l ] + = 

eigenvalues ±1. They are a linear combinations of the 
spinor operators T n |15j : 

2N 

xt (uj) = ^2yn(uj)T n , (6) 

n=l 

where functions y n (w) are known, e.g. for free boundaries 
see Ref . [5] . If x is the "ordering" direction in a sense that 
the Pauli matrix |15j measures the spin at site n, then 



In 



-1 — Pn-l&n Tin — Pn-l°~n 



(7) 
,N 



where P = 1 and P n = H" =1 (-crij), j = 1, 
and ctJj are Pauli operators |15j . In the customery lan- 
guage, P n is a string of disordering operators because 
(— <rj) reverses the spin at site j. We emphasise that 
this method does not give the X(ujj) vaccum | $) (the 
maximal eigenvector, since 7(0;) > 0) constructively, un- 
like the SML case [TS] . The maximum eigenvalue of V is 
A max (N) = exp(l/2) Y,f=i li^jW) and the free energy 
f(N) per unit length in the (1, 0) direction is now 



f(N) 



1 N 

-^ 7 K(7V)), 



(8) 



valid for /11/12 > 0. This is evidently quite unlike ((3]) in 
form. The function 7(0;) is as 7(0;) in but with K\ 
and if 2 interchanged. The ojj(N),j = 1, ...,N are the 
solutions of the discretisation equation imposed by the 
strip geometry 



exp2iNu> = exp i (6i(u>) + 62(0})) . 



(9) 



In this approach, the interpretation of the angles 5i and 
82 is clear; they are phase shifts for reflection of lattice 
fermions from the surfaces and ([8|) is a consistency con- 
dition requiered by the presence of two such bounding 
surfaces at a finite separation N 9, 18]. The functions 
S(uj),Sj(uj),j — 1,2 are given by the same expressions as 
S(oj),5j(uj),j — 1,2, but with K\ and K 2 interchanged. 

With free boundary, a study of the solutions of ((Sj> 
and of the non-triviality of the eigenvectors which re- 
late to them mandates that we take the solutions in 
(0,7r). In fact, there is one solution in each interval 
[ic(j-l)/N,irj/N],j = l,...,N provided K{ > K 2 ,(T> 
T c ). But for K\ < K 2 ,(T < T c ), there is a solution 
as above with j = 2, . . . ,N. If j = 1, there is such a 
solution but only if N < (<5*)'(0). If we take the pair 
-uj 1 (N),uj 1 (N), then as we pass from N < (S*)'(0) to 
N > (<$*)' (0), there is a bifurcation to a pair of imaginary 
solutions u) = ±iv(N), one of which is to be taken in the 
solution set. For v(N) > 0, up to order exp(— 4^7(0)) we 
have 



v(N) ~ 7(0) - 2 (sinh2if 2 sinh2ifi 



sinh 7 (0)e- 2A, ^W, 
(10) 



with 7(0) = 2{K\ - K$) (see Eq. (gj). The corresponding 
"j(iv(N)) is given up to order exp(— 27V7(0)) 



j(iv(N)) ~ 2smh2i4r^sinh7(0)e 



-AT^(O) 



(11) 



which clearly demonstrates asymptotic degeneracy as 
N —> 00 [9THTj. The bifurcation of a solution correspond- 
ing to j — 1 is a mathematical manifestation of the shift of 
a pseudocritical temperature in a film with respect to the 
bulk critical temperature T c . Recall that in a confined 
two-dimensional Ising system there is no true ordering 
phase transition; it becomes exponentially rounded and 
shifted towards lower temperatures. Such a shift, which 
is ~ (1/N), has been predicted from a scaling analysis by 
Nakanishi and Fisher [22] . 

In order to extract the Casimir free energy from (j6|), we 
use a contour integral method with a suitable summation 
kernel. The first use of a summation kernel in problems of 
this type, somewhat tangentially, is in Onsager's calcula- 
tion of surface tension [2Tj . There, the wavenumbers are 
equally spaced. Further, there are applications in stochas- 
tic processes |23 , and, moreover, the general problem was 
considered in [9], but in a different context. The contour 
integral is 



1 N 

-Y 

2 ^ 

3=1 



7MA0) = ^-_£ dtJ 7M 



1 dF, 



N 



Fn(u>) dui 



R (12) 



where R — — (l/2)(7(0) + 7(71")), with summation kernel 
F N (w) = e 2lNuj — e % ( s ( u »', the integration contour C sur- 
rounds every zero of F^(uj) in the strip — 7r < Reuj < n, 
with special attention to the behavior at u> = ±7r and 
U3 = 0. None of these is an allowed value, since the asso- 
ciated eigenfunction is trivial. The terms 7(0) and j(tt) 
do not report in the thermal Casimir force as they are 
N- independent. The development of the integral in (fl~2j) 
will be given elswhere. After extracting the bulk term, we 
find that the rather different representation are indeed the 
same. Since ((5]) has a scaling limit (by direct construc- 
tion), so does the Casimir free energy in the equivalent 
form 



fcas(N) 



1 . N . N 

E7WW))+ y 

.7=1 • y_7r 



7(w)dw (13) 



a result which would be difficult to extract from a direct 
assault on ([5]). We wish to assess the contribution of a 
single term in the sum to (| 13[) . Provided this has a scaling 
limit, then it is meaningful to compare it with ([5]). 

In particular, it is pertinent to compare the contribu- 
tion of (l/2)^y(iv(N)) coming from a complex mode to the 
Casimir force in the scaling limit as calculated numerically 
from © (see Eq. (3) in Ref. [12]). This is shown in Fig. [2] 
for K\ = K 2 and free boundaries. The Casimir force is the 
negative of derivative of the Casimir free energy with re- 
spect to N. Taking the scaling limit as N — > 00, 7(0) — > 0, 
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such that x = Nj(0)sgn(T — T c ) is fixed of the contribu- 
tion to the Casimir force from the imaginary mode, one 
obtains 



T c — — 

•> Cas — 



sinh2i^ 1 *(T c ) 1 



z 2 \fx r 



N 2 {\x\ 



z 2 )' 



for x < 0. The function z(x) is calculated from 
exp(— 2z) = (x + z)/{x — z). 



(14) 



(15) 



For x < — 1 there is one solution, for which the denomi- 
nator of (|T4| is strictly positive. For — 1 < x < ( and 
incidently for x > ) (fl5|) has no real solution. As one can 
see from Fig.[2]the imaginary wave number yields the dom- 
inant contribution to the scaling function Q(x) of the crit- 
ical Casimir force which is negative. Moreover, its magni- 
tude is larger than that of O. From Fig. [2 one can see that 
the surface state contribution governs the leading decay of 
the critical Casimir interaction; it can be shown that (fl"4"|) 
as x -> -oo is -{l/N 2 )(sinh2Kl(T c ))- 1 x 2 e~ x . At the 
bifurcation point x = —1, (1/2)7(1^1 ) goes to zero and is 
continuous in its argument. The contribution coming from 
the real mode uj\ is large and positive for — 1 < x < ( 
see blue dashed curve in the inset in Fig. [2]) and has to 
be compensated in order to render the negative critical 
Casimir force in this interval. This contribution is given 
by (TUT) with z replaced by iz, when now z is given by a 
solution of tan(z) = z/x for — 1 < x < 1. The Casimir 
force for (++) BCs may be obtained from the case of the 
free BCs by applying duality [16]. In this case there is 
an asymptotic degeneracy associated with a surface mode 
in the T > T c sector. Eqs. flT4l and (fTS"]) can be gener- 
alized to the case with the surface fields with hi and /12 
such that h\hi > [24]. This is an interesting matter 
because surface fields maybe adjusted to produce both at- 
tractive and repulsive Casimir forces. Equally, considering 
the case h\h2 < and scaling around the wetting transi- 
tion but not the bulk critical point, we can produced novel 
behavior. These are rather detail calculations and results 
which seemed to the authors to make separate publication 
advisable. 



We now give a more intuitive account of the asymp- 
totic degeneracy. For N large enough, but finite, we will 
have states in the strip which are either largely positively 
or largely negatively magnetized. At the transfer matrix 
level, called these states | ffi) and | 0). We can choose the 
phases so that the parity operator Pm | ©) =| 9) (The 
operator Pn is given by Eq. with n = N; note that 
Pjjf = 1. It is called the parity operator because it re- 
verses all x-quantised spins in the column simultaneously; 
it is an invariant of the Hamiltonian.) Then the fact that 
P N I $) =| $) implies that | $) = (| ©)+ | 9»/\/2 and 








-0.2 



Fig. 2: (color online) The surface mode contribution to the 
scaling function of the critical Casimir force Q c — N 2 ^^^ 
for the isotropic (Ki — K2) Ising strip with free boundary 
conditions (dashed green line) shown together with the total 
scaling function Q(x) (solid black line). The red dot-dash line 
is the result of the substraction of the part coming from the 
imaginary wave number k — iv(N) from the 0(x). The surface 
mode contribution (the dashed green line) is equal to at x — 
— 1. The dashed blue line shows the contribution to the scaling 
function of the critical Casimir force coming from the real wave 
number uii for — 1 < x < 0. 



further that X*(iv(N)) | $) = (| ©)- | ©)) /a/2, since 
P N X^(iv(N)) I $) = - I $). It follows that 



(1 + XHyu)) I $) 
V2 



e 



(l-xt(itt)) 1$) 

V2 



(16) 

The Euclidean "evolution" of these states is instructive: 
consider the matrix elements (© | V n \ ©) = (1 + 
e -n 7 (™W))/2 and (© I V n | ©) = (1 - e -nj^v(N))y 2 . 
Now ~f(iv(N)) - 7(0)exp(-A^7(0)) and if 717(0) > 
exp(7V7(0)), the transition matrix element tends to 1/2, 
that is, the system dephases on this length scale. Thus 
the dominant configurations correspond to successive pos- 
tively and negatively magnetized regions separated by do- 
main walls running across the strip as in the Fisher Priv- 
man theory (see Fig. lc in Ref. [12]). On the other hand, 
if n = 1, then (9 | V" | ©) ~ j(iv)/2 so that magnetic 
domain reversal is improbable but strictly not impossible 
if N < 00. Since (© | V n | ©) = (9 | V n \ ©), we see 



p-4 



Surface states and the thermal Casimir interaction 



that the wall interactions of the ® and phases are the 
same, an intuitively obvious result of symmetry. This is 
the physical reason for the smallness of j(iv(N)) in (TTT1) . 
In ([5]), with uji — iv(N) we have j/2n-i ~ -^o ex P(~ nv(N)) 
and y 2n ~ iN Q exp(-(N - n)v(N)) with v(N) ~ 7(0) 
(Nq is a suitable normalization constant). The imaginary 
mode qualifies as a surface state because of the exponen- 
tial attenuation away from the surface. This is because 
7(0) > 0, except at the critical point, which has a new at- 
tribute: it is also a de-pinning point for the surface state, 
simultaneously at both edges as 7(0) — > 0. In the scal- 
ing regime, the two exponentially decaying tails overlap, 
a physical reason to expect a significant contribution to 
the Casimir force. The bubble/domain wall, or interface 
Hamiltonian, idea has proved extremely useful in under- 
standing planar uniaxial models [25-27 . It, too, needs 
amplification to work with free boundaries: for modes 
with real ujj(N), path fluctuations are controlled by sur- 
face stiffness and by repulsion from the edges, as expected. 
For imaginary mode we have a random succession of do- 
main walls which cross the strip, pointing on average in 
the (0, 1) direction following Fisher-Privman. Within a 
stripe with positive magnetization, the state in a column is 
I ©) and it is essentially proved that the expectation value 
(© I r 2 „_i I ®) = t/2n-i- Now r 2 „_i consists of a string 
of disorder operators (the so-called the Jordan- Wigner tail 
see 0), thus it creats a column of misfit bonds with which 
a path may be associated. This path connects the end of 
the tail at the site n to any point on the boundary. Thus 
2/271-1 ~ exp(— nr(0, 1)) where r(0, 1) = 7(0) = v(N) is 
the surface tension in the (0, 1) direction. 

In summary, in this Letter we have obtained the follow- 
ing new results: (1.) Strong asymmetry of the Casimir 
force, as shown in Fig. [3J is related to the appearance in 
the spectrum of the transfer matrix of bound states, which 
are also surface states. (2.) Surface states produce the 
asymptotic degeneracy of the transfer matrix advocated 
as a diagnostic for phase transition by Kac [10]. They 
also play a key role in the realization of Fisher-Privman 
theory in these systems. The critical point has an addi- 
tional, new characterisation as an unbinding transition. 
We interpret the bound state at a microscopic level as a 
linear combination of block spin reversal operators. (3.) 
Our results require a significant modification of the inter- 
face Hamiltonian concept for these systems. 

Occurence of surface states depends on the boundary 
conditions of the strip. In a strip with free boundaries 
the asymptotic degeneracy occurs because the coexistence 
of two pseudophases in a finite system remains at the 
same line in the thermodynamic space (vanishing bulk 
field h = 0) as in the bulk system. Surface states form 
because free boundaries effectively attract domain walls 
due to the missing neighbours effect. This is analogous 
to the case of surfaces with wetting boundary conditions 
below T c , where the effective attraction due to the missing 
neighbours is amplified by the action of the surface fields. 
As a result, the thin film can be formed near the surface 



which then undergoes unbinding at the wetting transition. 
This argument carries over to three dimensional systems. 
For (+, +) boundary conditions an imaginary wavenumber 
and the associated asymptotic degeneracy occurs above 
T c ; there is no asymptotic degeneracy below T c because 
the surface breaks the symmetry and pseudo coexistence 
is shifted away from the h = line. For this boundary 
conditions, the imaginary mode gives rise to the critical 
adsorption which is characterized by the magnetization 
profile which is enhanced near both surfaces and decays 
to zero in the middle of the strip. 
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